DECAY OF POSITIVE WAVES FOR n x n HYPERBOLIC 
SYSTEMS OF BALANCE LAWS 



PAOLA GOATIN AND LAURENT GOSSE 

Abstract. We prove Olemik-type decay estimates for entropy solutions of 
n X n strictly hyperbolic systems of balance laws built out of a wave-front 
tracking procedure inside which the source term is treated as a nonconservative 
product localized on a discrete lattice. 



1. Introduction 

A classical result proved by Oleinik [18j for strictly convex scalar conservation 
laws in one space dimension shows that the density of positive waves decays in time 
like 0(l/t), see also [14]. More precisely, if we consider the scalar equation 

(1.1) ut + /(«)x = 0, u{t = 0,.)=u o eL°°(M), 

with f"(u) > k > 0, we have that every entropy-admissible solution satisfies 

(1.2) u(t,y) - u(t,x) <^—^ for alii > 0, x < y. 

and therefore has locally bounded variation (see also [S], Theorem 11.2.2). Con- 
versely, if u = u(t, x) is a weak solution satisfying (|1.2p . then u is entropy admissible. 

The same estimate as in (|1.2[) has been recovered for the Riemann coordinates 
of a particular 2x2 system [4] and for n x n genuinely nonlinear systems belonging 
to the Temple class 7 . However, one cannot expect such a result to remain valid 
for general n x n systems, even assuming all characteristic fields being genuinely 
nonlinear. Indeed, interactions among existing shocks may generate rarefactions as 
time increases. Decay estimates must therefore take into account for the generation 
of new positive waves due to interactions. Results in this direction were proved by 
Liu [17] in the case of approximate solutions constructed by Glimm's scheme, [5], 
and by Bressan and Colombo [5] for exact solutions obtained as limits of front 
tracking appoximations for n x n homogeneous systems, |19j . This in turn yields 
uniqueness of solutions satisfying the Oleinik entropy condition [SJ [7J [TO] . From the 
point of view of practical applications, such one-sided estimates are very useful for 
instance in the context of multiphase Geometric Optics computations, [12], or local 
error estimates, [20) and asymptotic behaviour of entropy solutions. 

In this paper we are interested in extending Olemik-type estimates on positive 
waves to quasilinear systems of balance laws. More precisely, we shall deal with the 
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Cauchy problem for the following n x n system of equations 
(1.3) u t + f(u) x = g(x,u), xeR,t>0, 



endowed with a (suitably small) initial data u a G L 1 n BV(JR; JR n ). Here u(t, x) G 
M n is the unknown function, / : O — » M n is a smooth C 2 vector field defined on 
an open neighborhood f2 of the origin in M n . We will assume that the system 
(|1.3p is strictly hyperbolic, with each characteristic field either genuinely nonlinear 
or linearly degenerate in the sense of Lax [T5] . Moreover we assume the following 
Caratheodory-type conditions for the source term g: 

(PI) g : M x fl — > M n is measurable w.r.t. x, for any u G f2, and is C 2 w.r.t. m, 
for any x G 5?; 

(P2) \\g(x, -)||c 2 is bounded over O, uniformly in a; G M; 

(P3) there exists a function uu G L 1 (J?) n L co (M) such that |g(x,u)| < u>(x), 
\\y u g(x,u)\\ < ui(x) for all (x, u) G K x SI. 
In addition, we require that a non-resonance condition holds, that is, the character- 
istic speeds are bounded away from zero: for some p G {1, . . . , n} and some c > 
one has 



for all u G fl, where Aj(w) denote the eigenvalues of the Jacobian matrix Df(u). 

Under these assumptions, it was proved in py that there exists a family of entropy 
weak solutions to (|1.3|) continuously depending on the initial data. More precisely, 
if the L 1 -norm of u is small enough, there exist a closed domain T> C L 1 (JR; M n ) 
of functions with sufficiently small total variation, a constant L and a unique semi- 
group P : [0, +oo) xD-tP with the properties: 

(i) For all u,v G V and t, s > one has \\P s u — P^Hl 1 < i(|t — + — 

(ii) For all u a G T> the function u(t, ■) = Ptu a is a weak entropy solution of the 
Cauchy problem (jl.3p . w(i = 0, .) = u . 

Under the above assumptions we aim at showing that, for genuinely nonlinear 
characteristic fields, an Olemik type estimate on the decay of positive waves holds, 
which takes into account not only new waves generated by interactions but even 
the contribution of the source term. A careful statement of these results requires 
some notations, [H El [6] . 

As usual, let A(u) = Df(u) be the Jacobian matrix of /, and call Aj(tt), k(u), 
ri(u) respectively the eigenvalues and the left and right eigenvctors of A(u). Let u : 
M — * Q have bounded variation in x and satisfy (|1.3j) with g = 0. The distributional 
derivative fj, = D x u is a vector measure. For i — 1, . . . , n we can now define fi l as 



where h(x) = li(u(x)) at points where u is continuous, while U(x a ) is some vector 
which satisfies 



where with a l a we denote the strength of the i-th wave generated by the resolution 
of the corresponding discontinuity in x a . 



(1.4) 






\k(x a ) - k(u(x a ))\ = 0(1) ■ \u(x a +) - u(x, 
k(x a ) ■ (u(x a +) - u(x a -)) = a l a , 



a 



)l 
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We denote by // the positive and negative parts of fj,, then we have 
I*' =!*'+-!**-, \» i \=( J , i+ +» i -. 
The total strength of waves in u is defined as 

n 

(1.5) V(u) =£><(«), Vi(u) = W\(M), 

i=l 

while the interaction potential is defined in terms of product measures on M 2 : 
(1-6) 

Q(«) = E (1^1 x lMi)({(^y) ^ < y}) + E 0**" x kl) ({(*,!/) = * ^ i/})> 

where (?A/" denotes the set of genuinely nonlinear families. 
Now we are ready to state our main result in the case g ^ 0: 

Theorem 1.1. (Decay of positive waves) Let the system (| 1.3$ be strictly hyper- 
bolic and let the i-th characteristic field be genuinely non-linear. Then there exists a 
constant C depending solely on f such that for every < s < t and every solution u 
with small total variation obtained as limit of wave-front tracking approximations, 
the measure fj, l t + of i-waves in u(t, •) satisfies 

(i.7) tf(J) < c ■ + Q(s ) _ Q( t) + v( Uo ) • |M| L1 

for every Borel set J C 1R. 

Of course, we tacitly assume (P1)-(P3) and (|1.4p throughout all the text. 

2. Wave-front tracking with zero-waves 

In this section we briefly recall the construction of wave-front tracking approxi- 
mations as stated in pQ. We start with the definition of the /i-Riemann solver. For 
small h > we introduce the map 



$ h (x ,u) = f 



-l 



h 

f(u) + I g(x a + s, u)ds 



(note that / is invertible due to (|1.4p ), that approximates the flow of the stationary 
equation associated to ()1.3|) . Consider now the Riemann problem with initial states 

ui if x < x , 
u r if x > x Q . 



(2.1) u(0,x) 



To locally render the source term's effects, a stationary discontinuity is introduced 
along the line x = x , that is, a wave whose speed is equal to zero; it will be referred 
to as a zero-wave. An ft,-Riemann solver for (| 1 . 3[) — (|2.1[) has been defined in [1] 
as a self-similar function u(t, x) = Rh{{x — x a )/t; ui, u r ) as follows: 

(a) there exist two states u~, u + which satisfy u + = $/j(a; ,u _ ); 

(b) u(t,x) coincides, on the set {t > 0,x < x Q }, with the solution to the 
homogeneous Riemann problem with initial values ui, u~ and, on the set 
{t > 0, x > x D }, with the solution to the homogeneous Riemann problem 
with initial values u + , u r ; 

(c) the Riemann problem between ui and u~ is solved only by waves with 
negative speed (i.e. of the families 1, . . . ,p); 
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(d) the Riemann problem between u + and u r is solved only by waves with 
positive speed (i.e. of the families p + 1, . . . , n). 

This clearly shares a lot of common features with the nonconservative Riemann 
problems studied in [16] . Let now e, h > be given: an e, /i-approximate solution 
of (II. 3j) is constructed as follows. First of all, the source term is localized by means 
of a Dirac comb along zero- waves located on the lattice x = jh, j 6 (— 4, 7^7) n 

(2.2) u t + /(u)x = ft ^ g(x, u).S(x - jh), 

j 

where <5 stands for the Dirac measure concentrated on x — 0. 

Given the initial data u Q , we deduce a piecewise constant approximation u(0, •) 
without increasing its BV-norm and u(t, x) is constructed, for small t, by applying 
the ft-Riemann solver at every point x = jh, and by solving the remaining discon- 
tinuities in u(0, •) using a classical homogeneous Riemann solver (rarefaction waves 
are discretized following [3 J: for a fixed small parameter v, each rarefaction of size 
a is divided, at its starting time, into m = [2] + 1 wave fronts of size tr/m < v). 

At every interaction point, a new Riemann problem arises. Notice that because 
of their null speed, zero-waves cannot interact among each other. In order to keep 
finite the total number of wave- fronts, two distinct procedures are used for solving a 
Riemann problem: an accurate method, which possibly creates several new fronts, 
and a simplified method , which minimizes the number of new wave-fronts. For a 
detailed description, as well as the proof of the stability of the algorithm, we refer 
the reader to [TJ [3] . 

The approximate solution can have four types of jumps: shocks (or contact 
discontinuities), rarefaction fronts, non-physical waves and zero- waves: J = S U 
1Z U MV U Z. A priori bounds on the functions u e ' h are obtained by modifying 
slightly the Glimm's functionals [S] in order to keep track of the zero-waves, 

(2-3) V(t)=J2W«\= E 1^1. 

aej aeSUKUAfVUZ 

(2.4) Q(u(t))= £ \<r a a p \<V{t) 3 , 

a,f3£A 

measuring respectively the total wave strengths and the interaction potential in 
u(t, ■). In particular, if a € 2 then the strength of the wave located in x a = j a h 
can be measured by means of 

f h 

(2.5) <J a — I uj{j a h + s)ds. 

Jo 

In (12. 4[) we have denoted by A an extended set of approaching waves. As usual we 
call k a the family of the front located at x a , with size o~ a . More precisely, a couple 
of wave- fronts of families k a , kp, located at x a < xp, belongs to A in any of the 
following cases: 

- if none of the two is a zero- wave, either k a > hp, or else k a = kp and at 
least one of them is a genuinely nonlinear shock, [9] ; 

- if a is a zero-wave and j3 is a physical one, kp < p; 

- if j3 is a zero-wave and a is a physical one, k a > P- 
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Notice that for some C > 1 there holds 

i[lM| L i (J) +Tot.Var.u(V)] < V(t) < c[||w|| L i (I) + Tot.Var.«(*, •) 
where the interval / is defined by 

i = U u* h > v° + w = U fo' ft > 0' + w ■ 

aez J'e(-&,&)n* 
Passing to the limit as e — > 0, /i > fixed, one has: 

(i) the total variation of u e ' h (t, •) remains uniformly bounded; 

(ii) the maximum size of rarefaction fronts approaches zero; 

(iii) the total strength of all non-physical waves approaches zero. 

By (i), Helly's theorem guarantees the existence of a subsequence strongly con- 
vergent in Lmoc. By (ii) and (iii), this limit provides a weak solution to (|2.2[) in 
agreement with nonconservative theories, [16 . At this stage, we can extract again 
a subsequence u hi which converges to some function u in L 1 ; oc and solves (II. 3|) in 
the usual weak distributional sense. 

In the sequel, we will need a semicontinuity property of Glimm's functionals. 
For h > fixed, let the total strength of waves and the interaction potential 
Qh be as in pQ, Section 4.1, that is, they are defined as in (|1.5[) - (|1.6p . but including 
the zero- waves. We have 

(2.6) V & (u)=V(u) + |M| L i, 

(2.7) Q(«) < Q h < Q(u) + |M| L1 ■ V(«). 

Proceeding as in [5] we recover the lower semicontinuity of the functionals Q/j and 
T 'h(u)=V \(u) + C Qh(u), C > 0, on a domain V of the form 

V=\u e L 1 n BV(fi; M n ), T h (u) < 7}, 7 small enough, 

(see [3], Theorem 10.1): 

Theorem 2.1. (Lower semicontinuity of the Glimm functionals) There 
exists a choice of the constants C ,7 > such that, ifY^u) = ~Vh{u) + C Clh(u) < 
7, then for any sequence of functions u u G D , u u — > u in L 1 as v — -> 00, one has 

(2.8) Q h (u) <IiminfQ fc (« w ), 

u — *oo 

(2.9) T h (u) <liminfT h (u„). 

u— *oo 

Moreover, for every finite union of open intervals J = I\ U . . . U Is there holds: 

(2.10) fi l± (J) + C Q h (u) <liminf (fj^ (J) + C Q h (u u )) , i=l,...,n. 

v — ^00 

3. Proof of Theorem 11.11 

I. By Lipschitz continuous dependence of the trajectories it is not restrictive 
to assume s = 0, t = T. We will consider a particular converging sequence u u ' h of 
e„, /i-approximate solutions with the following properties: 

(i) each rarefaction front x a travels with the characteristic speed of the state 
on the right: 

x a = >^k a (u(x a +)); 
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(ii) each shock x a travels with a speed strictly contained between the right and 
the left characteristic speeds: 

X ka (u(x a +)) < x a < X ka (u(x a -)); 

(hi) calling 7V„ the number of jumps in u v ,h = u u ' (0, ■), as v — > oo one has 

(3.1) e„ -> e u N v -> 0; 
(iv) the interaction potential satisfies 

(3.2) Q h {u v ' h (0, •)) -> Qh(u ) as v -> oo. 

Such a sequence can be constructed as explained in 3„ proof of Lemma 10.2 (p. 205). 

Let u = u(t,x) be a piecewise constant e, h- approximate solution constructed 
via front-tracking approximation (we shall drop the e, h superscripts since there is 
no ambiguity). As usual, by (generalized) i- characteristic we mean an absolutely 
continuous curve x = x(t) such that, [8], 

x(t) 6 [\i{u(t, x+)), Xi(u(t, x-))] a.e. t > 

By £ i—* y I (t; x) we denote the minimal z-characteristic passing through x, at time T. 
Because of (|1.4|l the presence of zero-waves does not affect the usual construction. 
Let now I = [a, b[ be any half-open interval, and define 

I{t) = [y\t;a), y \t;b)[= [«(*), &(*)[■ 

We seek an estimate of the amount of positive i-waves in the the approximate 
solution u(T, •) contained in I. We recall that k a stands for the family of the front 
located at x a , with size a a . For a genuinely nonlinear family, the size of the jump 
can be measured like 

(3.3) a a = X ka (u(x a +)) - X ka (u(x a -)), 
while the size of a zero- wave is still given by (|2.5[) . Define 

m(t) = b(t) - a(t). 

By (I3.3[) and the Lipschitz continuity of the map u Xi(u) we deduce that 
m(t) = Xi(u(t,b(t))) -Xi(u(t,a(t))) 

(3.4) = M(t) + 0(l)(e + K(t)) 

for a.e. t. The Landau symbol stands for a quantity whose modulus is uniformly 
bounded. We see that 

M(t)= J2 * a =i4(m) 

is the total amount of (signed) i-waves in u(t, •) contained in I{t), while 

K(t)= E Ki = EKl( 7 (*)) + / ^- dx 

stands for the total strength of waves of families ^ i inside I(t), zero- waves included. 
To estimate the contribution of the term Kit) in (|3.4|) we introduce 

*(t) = E 0k a (t,x a (t)) ■ \a a \ < V h («(t)), 

k a ^i 
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where 

if x < a{t) 

<•,(/..<> = <{ ~7T?~r i£xe[a{t),b(t)[ 
if x > b{t) 

if x < a(i) 

",(/•<■)=<{ -^f- if* e[a(t),&(t)[ 
if ar > 

in the cases j < i or j > i respectively. Roughly speaking, $(t) represents the 
cumulated strength of the waves which do not approach the interval I(t). By strict 
hyperbolicity, we can expect it to grow with time. Observe that $ is piecewise 
Lipschitz continuous with a finite number of discontinuities occurring at interaction 
times, where it may decrease at most of: 

(3.5) 4>(t+) - *(r-) = 0(l)[Q(u(r-)) - Q(u(r+))}. 
Following [3], we assume that there holds for some c > 0, 

(3.6) \Xi(u) — Xi(v)\ < c , \\i(u) - Xk a (v)\ > 2c , 

for every couple of states u,v and every k a ^ i. Outside interaction times $ is 
non-decreasing; indeed, we have 



k a <i,x a el(t) 



2^ \ a ° 



b — x a {b~x a )m 



{x a — a)rh 



k a >i,xa,el(t) 



in 



> V ■ 



^>-"> m(t) 



thanks to the system's strict hyperbolicity and the non-resonance condition l|1.4[) . 
In particular, observe that c Q < c in ()3.6|) . The above estimate yields the bound 
valid for all but finitely many times t: 

(3.7) K(t) < —<f>(t)m(t). 

Co 

We notice again, as in [TJ [^] , that there's a need for a completely different theory 
in order to tackle resonant cases, see also |13) . 

Concerning the term M(t), observe that it can change only when an interaction 
occurs within the interval [a(t),b(t)]. In this case, one has 

M(r+) - M(r-) = 0(l)[Q(u(r-)) - Q(u(r+))}. 

This yields an estimate of the form 

(3.8) M(T) - M(t) = 0{l) [Q(«(r-)) - Q(«(r+))] 

rer 
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where the summation extends over all times r G]0, T] where an interaction occurs 
inside [a(r),6(r)]. Inserting the estimates (|3.7p — (|3.8j) in (13. 4|) we obtain 

(3.9) m(f) + C*(t)m(t) > M(T) - C + ^ |AQ(r)|J , 

for some constant C and a.e. t. We now observe that m is a continuous, piecewise 
linear function of t, and $ is uniformly bounded. It can decrease only at interaction 
times, where (|3.5[) holds. Hence its total variation in uniformly bounded and for 
some constant K Q we have the estimate 

(3.10) \ §(t)dt < K Q . 

Jo 

As in [3], from (|3.9p we deduce the decay estimate: 

(3.11) M(T) < 2e CK ° (~~jT~~^j + 2Ce + 2C l A ^( T )l- 

II. Repeating the above process for any finite number S of disjoint half-open 
intervals I s = [a s ,b s [ we obtain 

(3.12) J2 M *( T ) ^ C ' ( E — f^ 1 + Se + [Q(«(0)) - Q(«(T))] j , 

s=l Vs=l / 

for some constant C" independent of and of the particular e, /i-approximate solu- 
tion. Here we have used the notation 

M S (T)= Yl a « = A*t([o..6-D 

to denote the sum of (signed) strength of all i-waves in u(T, •) contained in the 
interval [a s , b s [. 

III. Let us consider now any open interval ]a, b[. Let N be the number of 
i-shocks of the first generation in the front-tracking e, /i-approximate solution u, 
as defined in [3], Chapter 7. We can thus construct half-open intervals I s = [a s ,b s [, 
s = 1, . . . , S < N + 1, such that the following holds (see [3J, Chapter 10): 

• Every i-rarefaction front in u(T, ■) contained in ]a, b[ falls inside one of the 
intervals / s . 

• No i-shock front of the first generation fall inside any of the intervals I s . 
Calling /4+ the measure of positive i-waves in u(T, •), the above properties imply 

(3.13) /4+(K b[) = Y, W) + 0(1) • [Qfc(«(0)) - Qa(u(T))] + 0(1) • r(e) 

where Q/j is the interaction potential introduced in [T], Section 4.1. 

Indeed the only negative i-waves contained in (J 7 S must have generation order 
> 2, originating from interactions during the time interval ]0, T]. The total strength 
of these negative i-waves is bounded by the decrease in the interaction potential Q. 
The last term on the right hand side of ()3 . 1 3|) tends to zero as e does and comes 
from the difference between Q(u) and Qfc(u): in the latter there are no non-physical 
fronts and all the countable ft,-Riemann problems are solved (whereas before only 
zero- waves inside the interval (— 1/e, 1/e) were considered). 
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Altogether, lETT2"jl and ([3713]) yield 

(3.14) /4+(]a, b[) < C" + (N + l)r(e) + [Q h (u(0)) - Q h (u(T))} 

for some constant C" independent of e. 

IV. For h > fixed, we now consider a sequence of e u , /i-approximate solutions 
satisfying the properties (i) — (iv) stated at the beginning of this section. It is 
clearly not restrictive to take C" > C in (|3. 14[> . where C is the (big) constant in 
Theorem O Using flU]), IpTO]) . (j3~TIj) . (31} and we obtain 

Mt + (KM) 

< liminf (mL + t(]MD + C"Q h (u v (T))) - C"Q h (u(T)) 

< C" • liminf ( h —^ + (i\T„ + l)r(e) + Q h K(0)) ) - C"Q h (u(T)) 
QfcKO)) - Q h (u(T)) 



T 
T 



Q( U (0))-QKT)) + [[o;|| Ll -VK0)) 



Since the last term is independent of h this proves (| 1 . T|) in case J is an open interval. 
The same arguments can be used in the case where J is a finite collection of open 
intervals. Since fi 1 is a bounded Radon measure, the estimate fll.7p holds for every 
Borel set J and we are done. □ 

4. A SHORT COMMENT 

All these computations heavily rely on the restrictive assumption that the source 
term is dominated by a function uj G L 1 (J?) (the L°°(iR) bound being just a 
consequence of the smoothness of g). In the case of a convex scalar law, [IT], 

u t + f(u) x = g(u), u e L 1 nBV(2R), 



the stability of the approximation procedure for (|1.3[) can be obtained requiring 
only u> £ L°° (JR) . Hence one can follow the same canvas and (|3.4p simplifies a lot 
thanks to the obvious bound K(t) < \\u)\\x 1 °°m(t). Thus (|3.9|) boils down to 

m(t) + |M| L oom(i) > M(T), 

which leads in sharp contrast to an exponential bound of the type (fj, 1 — u x ): 

(4.1) f i 1 t + (J)<Ce^°° t - meaS ( J \ t>0. 

This agrees of course with simple computations, since in this context Oleinik's 
estimates can be derived from the Riccati differential equation: 

z t - g'(u)z + f"(u)z 2 = 0, z(0) = sup ( max (0, (u ) x ) X e H. 
Therefore one recovers f|4. 1 [) in the case C = 1/k and uj = Lip(g) since we have: 

„Lip(g).t Lip(g).t 

z(t) = —. ^< , t>0. 

' 1 ■ • Kb 



1 \~( ef^P^J-lA 
7m +K \ L lp (g) J 



MP) 

All in all, the L 1 (K) bound on u> expresses somehow the fact that the source 
has negligible effects outside a compact interval in M as pointed out in [5] p. 329. 
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Hence, by the strict hyperbolicity ensured by ()1.4p - (|3.6j) . waves exit this region 
after some time and then are ruled only by the convective process. 
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